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G.Breit's original paper of 1929 postulates the Breit equation as a correction to an earlier defective 
equation due to Eddington and Gaunt, containing a form of interaction suggested by Heisenberg and 
Pauli. We observe that manifestly covariant electromagnetic Two-Body Dirac equations previously 
obtained by us in the framework of Relativistic Constraint Mechanics reproduce the spectral results 
of the Breit equation but through an interaction structure that contains that of Eddington and 
Gaunt. By repeating for our equation the analysis that Breit used to demonstrate the superiority of 
his equation to that of Eddington and Gaunt, we show that the historically unfamiliar interaction 
structures of Two-Body Dirac equations (in Breit-like form) are just what is needed to correct the 
covariant Eddington Gaunt equation without resorting to Breit 's version of retardation. 
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I. INTRODUCTION 



Three score and seven years ago, Gregory Breit extended Dirac's spin-1/2 wave equation to a system of two charged 
particles [1] . He formed his equation by summing two free-particle Dirac Hamiltonians with an interaction obtained 
by substituting Dirac a's for velocities in the semi-relativistic clcctrodynamic interaction of Darwin:. 



= |«i • Pi + Pim-i +a 2 -p2+ [3 2 m 2 - ^[1 - • d 2 + cti ■ f a 2 ■ r)]| (1.1) 



Although successful for Breit's purpose - perturbative calculation of the electromagnetic bound-state spectrum 
of multi-electron atoms, Breit's equation turned out to be neither relativistic nor a well-defined wave equation. 
Nonetheless, ever since, people have applied Breit's equation in nuclear and particle physics in situations far from its 
electrodynamic origin. 

In the last 20 years, some of us [2-8] have tried to remedy this by constructing new fully covariant multiparticle 
relativistic quantum descriptions. L. Horwitz has had a hand in this effort as participant and critic. As he and F. 
Rohrlich pointed out in 1981 [9], for such descriptions, the two-body problem is essentially simpler than the generic n- 
body problem - reducing (just as in non-relativistic physics) to an effective one-body problem (most simply described 
in terms of effective variables introduced earlier by I.T.Todorov [10]). For most of these 20 years the authors of this 
paper have concentrated their efforts on using the special two-body equations to find sensible versions of them for 
interactions actually found in nature as opposed to "toys" like relativistic rotors and harmonic oscillators. We have 
applied them (with surprising success) to the phenomenological calculation of the q-q meson bound state spectrum 
[11,12]. 

But, of more fundamental but related importance, we have demonstrated that our equations reproduce the spectral 
consequences of the original work of Breit but in a surprising way: whereas Breit's equation yielded these results 
only through first-order perturbation theory, our "Two-Body Dirac equations" yield the same results through non- 
perturbative solution of a fully-relativistic quantum wave-equation [13,14]. For two relativistic spin-one-half particles 
interacting electromagnetically, our equations are given as a pair of compatible Dirac equations on a single 16- 
component wave-function: 

Drf = (tti • 7i + mi)V = (1.2) 
D 2 ip = (7T 2 • 72 + m 2 )ip = , (1.3) 

in which 

iTi = pi - Ai , i = 1,2 

and 

Ai = [l- \{G + G~ 1 )]pi + \{G- G~ 1 )p 2 - l -{dG- 72)72 

A 2 = [1 - i(G + G- X )\ V2 + \{G- G- r ) Pl +'-{dG- 71)71 

where G 2 = (1 — 2A/ w)^ 1 . w is the total cm. energy, while the world scalar potential A is a function of the covariant 
spacelike particle separation 

x" ± =x^ + P^{P-x) (1.4) 

which is perpendicular to the total four-momentum, P. For 0(a 4 ) electrodynamics, A = A(x±) = —a/r in which 
r = y f x'j_. In these coupled equations, the subscript i = 1,2 stands for the ith particle so that mi, and m 2 are the 
masses of the interacting fermions. The potentials A^ in this form of the Two-Body Dirac equations appear through 
minimal interaction substitutions on the particle four-momenta. They introduce mutual interactions as though each 
particle were in an external potential produced by the other. (Hence, we refer to these forms of the Two-Body Dirac 
equations as the "external potential forms" or the "minimal interaction forms".) The specific forms of the covariant 
spin-dependent terms in the interactions as well as the dependence of A on x± are consequences of the necessary 
compatibility of the two Dirac equations 

[£>i,£> 2 ]^ = 0. (1.5) 

(In most of our previous works, we have set out these equations and their compatibility in terms of S operators that 
result from the D's through multiplication by 754 's.) 
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Because we were originally able to obtain their versions for world scalar interaction by rigorously " taking the 
operator square roots" of two Klein-Gordon equations, to us they seem the natural two-body extensions of Dirac's 
one-body equation hence deserving the name "Two-Body Dirac equations". For the electromagnetic interaction of 
Eqs.(1.2-1.3), these equations possess a family of exact solutions for para-positronium with correct spectrum through 
0(a 4 )[13]. Moreover, they have not only been derived from quantum field theory (through a relativistic Lippmann- 
Schwingcr equation) as the " quantum- mechanical transform" of the Bcthc-Salpctcr equation [15, 16], but their spin- 
independent structures have been obtained as well from classical relativistic field theory through Whcclcr-Feynman 
electrodynamics. Most recently, Jallouli and Sazdjian have completed work begun for spinless charged particles by 
Todorov [10] more than 20 years ago, by obtaining our spin-one-half electrodynamic equations from the quantum- 
ficld-theoretic eikonal approximation [17]. 

For all their successes, people have been slow to adopt the Two-Body Dirac equations, preferring to stick with 
methods that on the one hand artificially truncate the fully-relativistic Bethe-Salpeter equation (and thereby produce 
well-known pathologies) or on the other abuse the Brcit equation by modifying its interaction structure in hopes that 
its nonperturbative problems will heal themselves. The perturbative success of Breit's equation in electrodynamics 
becomes a psychological barrier to use of the fully-relativistic nonperturbative equation. 



II. HISTORY 



The historical origins of this psychological effect are intriguing. To investigate them, we begin by noting that 
because the Two-Body Dirac equations are a compatible pair, they can be rearranged in a large number of equivalent 
forms. Among these is a "Breit-form" which in the cm. rest-frame reduces to the sum of two free Dirac Hamiltonian 
operators plus a covariant interaction. We set out this Breit Form in a previous volume of Foundations of Physics (in 
honor of F. Rohrlich)[12]. (The special properties that allow this equation to function as a non-singular version of the 
Breit equation were recently investigated by one of us in [18].) However, we have come to realize that the particular 
form of this interaction is part of an interesting controversy in the history of two-body equations. 

We decided to find out just what prompted Breit to formulate his famous equation. Reading his paper of 1929, 
we saw that his work was actually proposed as the correction of an older "defective" equation called by him "the 
Eddington-Gaunt equation" [19] 

= [ai ■p 1 +(3 1 m 1 +a 2 ■ p 2 + f3 2 m 2 - -(1 - a 1 • a 2 )]*, (2-1) 

r 

with an interaction structure first suggested by Heisenberg and Pauli (according to Brcit). We were startled to see 
that this equation was one already familiar to us as the lowest-order approximation to ours (referred to by us in [12] as 
"the familiar form for four- vector interactions without the Darwin piece"). In fact , in "Breit-like" form, our equation 
reads [12] 

= [d\ ■ pi + Pimi +a 2 -p 2 + [3 2 m 2 + w(l - exp[-<?(a;_i_)(l - «i • a 2 )])] , J / , (2.2) 

in which Q = ln[G] . 

Since our covariant equation automatically has correct retardation, apparently we have discovered in Eq.(2.2) a 
(higher order) corrected form of the (Heisenberg-Pauli) "Eddington-Gaunt" equation that has been "fixed" by a 
mechanism different from Breit's! 

What's going on here ? In fact closer reading of Breit and Eddington reveals a fundamental reason for all this. 
Breit judged the Eddington-Gaunt equation defective because of its failure after rearrangement as an equation for the 
"large components" of the wave function to reproduce the spinless Darwin interaction: 

a 1 

[i- - -(vi ■ V2 + V! ■ f v 2 ■ r)] (2.3) 

v £ 

and its consequent failure to yield a spectrum for He in agreement with experiment. But we recall that 
the Darwin interaction [20] is straightforwardly obtained by breaking the manifest covariance of J J J\GJ 2 = 
J dri J dT 2 XiX 2ll G[(xi — x 2 ) 2 ] in Wheeler-Feynman electrodynamics [21] by expanding the (1/2 advanced +1/2 
retarded) Green's function about the instantaneous limit. An integration by parts in the time leads to 

- -[(1 - V! ■ v 2 ) + \vx • (1 - f f) ■ v 2 ] (2.4) 
r 2 

with the 1/2 in the second "retardative term" coming from a Taylor expansion of the Green's function [22]. By redoing 
the Darwin interaction in terms of the cfs, Breit was able to restore the missing terms. On the other hand, if we read 
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Eddington's notorious paper [19] (and ignore his musings about the value of the fine-structure constant a), we find that 
he arrived at the — — (1 — d\ • d 2 ) interaction by performing the covariant four- velocity goes to (l,d?) substitution on 
the manifestly covariant JGJ coupling (that later became the heart of Wheeler- Feynman electrodynamics) in the form 
of ij 1 X2fj, f. In fact, this covariant form appears explicitly in our Eq.(2.2) where it is inherited (but correctly) directly 
from classical field theory (if we start from the Wheeler- Feynman approach [23]) or from quantum field theory in the 
form of the Bethe-Salpeter kernel in covariant Feynman Gauge [14]. This structure is a fundamental feature of the 
manifestly covariant approach! By comparing Eddington and Gaunt's Eq.(2.1) with our Eq.(2.2), we see immediately 
that where they went wrong was to stop at lowest order, thereby missing an all-important retardative recoil term! 

In detail, in Eq.(2.2), if we use the 7 matrix algebra to compute the exponential of matrix form to all orders, we 
find that 



w 



(l-exp[-g(x ± )(l-a 1 -a 2 )}) (2.5) 



= w eXp ^ G] [l(3cfe[g] + ch[3G}) + 751752 (3sh[G] - sh[3g}) 
+a 1 ■ a 2 {sh[Q] + sh[3G}) + <7 1 ■ a 2 (ch[g] - ch[30\)], 

a corrected form of Eq.(69) in [12]. If we eliminate g in terms of the potential A = —f, we find the striking result 
that to all orders in the potential: 

iu(l-exp[-0(aM_)(l-5i -a 2 )]) (2.6) 

A 2 A 3 1 

= .4(1 - d\ ■ a 2 ) (1 - u\ ■ S 2 ) - (1 - 751752 + a 2 - 5\ ■ a 2 ) 



w v " ' 1 w 2 (1 — 2 — ) 

From this, we immediately see that the perturbative dynamics (through 0(a 4 )) will be given by 

A 2 

A(l-ax-a 2 ) (I-0W2). (2.7) 

w 

Apparently, as a lowest-order perturbation, the second term of this must be dynamically equivalent to the term 

-ai ■ (1 - r r) ■ a 2 

r 2 



which when added to 



makes Breit's famous 



--(1 - di ■ d 2 ) 

r 



a 1 _ _ ^ _ 
- — [1 - -(ai • a 2 + ai ■ r a 2 ■ r)\ . 
r 2 



III. EQUIVALENCE? 

Not only must these terms be dynamically equivalent, but somehow our equation must have restored the content 
of the Darwin interaction that Breit discovered was spoiled by the Eddington-Gaunt equation. Breit uncovered this 
defect by performing a perturbative (in orders of 1/c) reduction to an equation on the upper- upper component of 
the 16-component wave function for both his equation and the Eddington-Gaunt equation. This rearranges the Breit 
equation as a set of corrections to the non-relativistic Schrodingcr equation with Coulomb potential. Expressed in 
the cm. frame this rearrangement becomes 



HiP = wip (3.1) 



in which w is the total cm. energy and 



v 2 (p 2 ) 2 v 2 (n 2 ) 
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p 2 1 1 

-a([(l ) p ■ (1 - ff) ■ p\ordered 

toito 2 r 2m2rri2r 



1,1 1 s , N 1 L 1 2 x ^ ,1 2 

2 m{ 4r mim 2 mim 2 



1 , 8?r ^ _ (?i • ct 2 3ai • f<T 2 • f v 

+ I (—5-0-1 ' o 2 5(r) + 3 5 ) . (3.2) 

4mim 2 3 r 6 r° 7 

To see what new features are contained in our truncated interaction Eq.(2.7), we perform the same type of semirel- 
ativistic reduction on our equation (see Appendix). We obtain a seemingly different Hamiltonian: 

P 2 (p 2 ) 2 P 2 (p 2 ) 

H = (m 1 + ^-- f4) + (m 2 + - ^i)+ 
2mi 8mJ 2m 2 8m 2 

! -{p 2 , 1/r} - — — ra/r 2 + -^—5{P) 



2toito 2 ' 2(mi + to 2 ) mirn 2 



r l , 1 1 .,,ir-p 



1L .. 1 2 _ .1 2 w . 

-7^T ' — 2 + )o"i + — + o 2 

4 r 6 ra\ mim 2 mim 2 

1 , 8tt ^ ^ cti • a 2 3cti • f<r 2 • f , 
+ i ( — 5-01 • cr 2 5(r) H 3 = ) . 3.3 

4mim 2 3 r 3 r 5 ' 

But, note that the r ■ p term and 6(f) term in the third line of Eq.(3.3) give equivalent expectation values (both 
contributing only to the ground state). Together they are equivalent to the first term of the third line of Eq.(3.2). 

Now, we come to the feature that led Brcit to discard the defective Eddington Gaunt equation in favor of his new 
equation: reproduction of the correct Darwin interaction. Here this issue boils down to the equivalence (or lack of 
equivalence) of the second lines of Eq.(3.3) and Eq.(3.2). As luck would have it, Schwingcr [24] has shown us just 
how to do this. In order to use the virial theorem to evaluate expectation values in his treatment of the positronium 
spectrum, Schwinger introduced a canonical transformation that turns the expectation value of the second line of the 
reduced Breit equation Eq.(3.2) into precisely the expectation value of the set of terms appearing as the second line in 
the reduced version Eq.(3.3) of our equation. In fact, we found some years ago that this transformation may be used to 
derive the relativistic Todorov equation of electrodynamics from classical electromagnetic field theory [25]. This takes 
care of the important spin-independent Darwin interaction that served as Breit's criterion for rejecting the Eddington 
Gaunt equation. But, what happened to the intriguing spin-structure of the new term? If we look into the details of 
our approximate reduction of our 16-component equation to a single 4-component equation on the upper-upper wave 
function (outlined in Appendix A), we find that the spin-spin term in the equation for the upper- upper component 
is exactly cancelled by a spin-spin term resulting from its coupling to the lower-lower component. Cancellations and 
simplifications of other spin dependences then yield the Fermi-Brcit form of our equation - Eq.(3.3). This establishes 
that the unfamiliar interaction of our Eq.(2.7) is perturbatively quantum-mechanically equivalent to that of Breit 
thereby providing a corrected form of the historical Eddington Gaunt equation. The lesson here is that manifestly 
covariant techniques yield manifestly covariant equations that produce correct spectral results if they contain the 
correct interactions. We have a situation in which some of the original investigators wrote down an interaction that 
was covariant (including Eddington's realization that the interparticle separation r should be something like \J 'x'j_ 
[19]) but incomplete, while another (Breit) replaced it with a semi-relativistic approximation to an interaction that 
was correct. 
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Appendix A - Semirelativistic reduction of Eq.(2.2) with approximate interaction Eq.(2.7) 

We write the sixteen component Dirac spinor as 



.1p4. 



(Al) 



in which the ipi are four component spinors. All the matrices that operate on this spinor are sixteen by sixteen. In 
the standard Dirac representation (the subscripts on the identity 1 give the dimensionality of the unit matrix) 



/Is \ (0 1s 



Pl751 = Pi 



1 8 

-Is 



(A2) 



(A3) 



75 \ (OU 

1 o 75 J' 75 = (i 4 o 



(A4) 



u 

-u 



(A5) 



/3l75l752 



75 
-75 



(A6) 



/?2752 7 51 = 



Op 
pO 



(A7) 



Using di = -f5i<Ti, our truncated Breit equation becomes 



a a 

[w - mi - m 2 H 1 - <?i • CT2)]V'i 



+P • <7 2 ^2 - P- <?l?p3 <?1 • 0*2^4 = 

r 



(A8) 



r a a /-, - m , 

w - mi + m 2 H 1 5-(l - (Ti • cr 2 ) V>2 



a 

+P ■ 0-2^1 - P ' 01^4 CTl ' 02^3 = 

r 



(A9) 



[to + mi - m 2 H 1 ^-(1 - <?i • ct 2 )]V'3 

r wr A 



_ _ _, _ a _ _ : 

+P ■ 02^4 - V ■ <J\W\ Cl ' 0"2"02 = 

r 



(A10) 



r a a /-, m / 

[to + mi + m 2 H 1 k(i - <J\ ■ (j 2 )\Wa 



_ _ . _ _ a _ _ . 

+p • (72^3 ~ P ' 01^2 0"1 ' 02^1 = 

r 



(All) 



G 



By using Coulomb variables in which p ~ a and 1/r ~ a we can obtain an expansion involving just the upper-upper 
component wave function ipi. The expansion we desire is one through order a 4 . To obtain this requires successive 
substitutions into Eq.(A.8) of expressions for the lower component wave functions given in Eqs.(A.9-ll) through 
appropriate orders. From Eq.(A.ll) we obtain 

ip 4 = — ; - ■ —[(a/r)ai -u 2 -p- S 2 ipz+p- aiip 2 ] (-4-12) 

w + mi + m 2 + a/r 

Note that we have ignored the a 2 spin-dependent term in the denominator since in the substitution into Eq.(A.8) it 
would have produced terms that are higher order than a 4 . We solve Eqs.(A.9-10) for ip 2 and 03 to obtain 

tt>2 = 7 TT7 TT 

(w + mi — m 2 + a/r)(w — mi + m 2 + a/r) 

x[(w + mi - m 2 + a/r)(p- c?iV>4 — p- a 2 ipi) + (a/r)Si ■ a 2 (p ■ aiipi - p- a 2 4>4,)] (A13) 

1 



tp3 



(w + mi — m 2 + a/r)(w — mi + m 2 + a/r) 



x[(w -mi+m 2 + a/r)(p- aiipi —p - a 2 ipt) + (a/r)5i ■ a 2 {p ■ aiipd. - p- <r 2 0i)] (^-14) 

In Eq.(A.8) we need 04 through order a 2 and thus in Eq.(A.12) we need ip 2 and 03 only through order a. This in 
turn implies that we can drop the 04 and interaction terms in Eqs.(A.13-14). Performing this substitution yields the 
following expression for 04 of appropriate order for use in Eq.(A.8): 

/ 1 \ a - - 1711 + m2 - - - - 1 / I A 1 K\ 

ip 4 = — ■ -o-i • <r 2 - — p ■ oip ■ a 2 \ipi (A.15) 

2(mi + m 2 ) r 2mim 2 

For the direct ^2 and 03 terms in Eq.(A.8) we need Eq.(A.13) and Eq.(A.14) through order a 3 , which in turn requires 
Eq.(A.15). This leads to 

, 1 P-<?1 ,0. _ _ ^ ^ ^ mi+m 2 

■02 = ■ ■ -[p-a 2 - — ■ ( &i -a 2 - p ■ oip ■ <J 2 — )\tpi 

w — mi+m 2 + a/r 2(mi + m 2 ) r 2mim 2 

(a/r) _ 
Amim 2 

, 1 r ^ - p-d 2 a _ _ _ _ _ _ mi + m 2 

V>3 = — ; 1 rP'ffi - TT, ! r(— o"! -a 2 - p ■ o\p ■ a 2 — )\ipi 

w + mi—m 2 + a/r 2{mi + m 2 ) r 2mim 2 

(n/r) a^^p-a^i. (A17) 



4mim2 

Combining all terms of Eq.(A.8) and replacing w by mi + m 2 in the a terms we obtain 



a a 2 



wipi = [mi + m 2 — ■ -(1 - cti • a 2 )]ipi 

r r z (m\ + m 2 ) 

. -. -. 1 r ^ -. f ot _ „ _ _ „ _ mi + m 2 

+P ■ o 2 ■ ■ r \p-a 2 - —. ■ r(-oi ■ <t 2 - p ■ oip ■ u 2 — )\ipi 

w — mi+m 2 + a/r 2{mi + m 2 ) r 2mim 2 

. -, -, 1 r - - P - 3 2 ,a _ _ _ _ _ _ mi + m 2 

+P ■ a\ ; ; -rXP ' <*\ - TT, ; ■ cr 2 - p ■ a x p ■ a 2 — ) Vi 

w + mi — m 2 + a/r 2(mi + m 2 ) r 2mim 2 

„ „ (a/r) _ _ _ _ _ _ (a/r) _ _ _ _ 

~P ■ °\ ■ &2P ■ 0"! "01 - p ■ (71- 1 • a 2 p ■ (72-01 

Amim 2 Amim 2 
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a _ _ a _ 1 , . 

H — o\ ■ (72 -7 ; T -0 '! ' CT 2 - -. P ■ <JiP ■ cr 2 Vi (A18) 

r l\m\ + m2)T Amirri2 

We omit most of the steps of the remaining reduction, involving powers of Pauli matrices, commenting on just two 
portions of the details. The first is that the spin-spin term in the first line cancels with the corresponding spin-spin 
term in the last line that results from the identity (p\ ■ (T2) 2 = 3 — 2<ti ■ 02- The second is that the semirelativistic p 1 
kinetic corrections not only result from terms at the end of the second and third lines that involve four p ■ <Ji factors 
but also from the terms at the beginning of each of those lines, in particular 

(p-a 2 ^ ■ —p-a 2 +p-vi — : : r p-Si)^. (A19) 

w — mi + m 2 + a/r w + mi — m 2 + a/r 

One brings the denominator through and operates on ipi using the lower order portions of wave equation Eq.(A.18). 
When we perform this operation and various other simplifications we obtain the semirelativistic reduction given in 
Eq.(3.3) in the text of the paper. 
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